THE SERIES LIMIT OF ]T fc l/[fclogfc(loglogfc) 2 ] 
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^v^j ' Abstract. The slowly converging series J3^L 3 l/[fclogfc(loglogfc) a ] is evalu- 

ated to 38.4067680928 at a = 2. After some initial terms, the infinite tail 
i of the sum is replaced by the integral of the associated interpolating function, 

which is available in simple analytic form. Biases that originate from the dif- 
ference between the smooth area under the function and the corresponding 
Riemann sum are corrected by standard means. The cases a = 3 and a = 4 
are computed in the same manner. 
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1. Aim and Scope 
We aim at a precise numerical evaluation of the series limit of 

1 



(i) c s £ 



, 3 Hogfc(loglogA:) 2 ' 



CN 

■ which has been estimated at C « 38.43 in the CRC tables [91 p. 42]. The presence 

0^ , of the square of the double logarithm in the denominator is just sufficient to achieve 

convergence; direct numerical summation is a futile strategy to estimate the series 
limit. This work addresses how the family of closely related 



■ Definition 1. 
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GN . (o) c [a) EE V 1 

O v ; ^ /clog fc(log log fc) 



is calculated by other — yet standard — methods of numerical analysis for a = 2 
to a = 4. 



2. Numerical strategy 

2.1. Romberg Integration. Some initial terms up to k — N are summed directly. 
Fig. [T] demonstrates the methodology for all larger indices: Starting at k = N + 
1, this could be made precise by adding the areas of the rectangular boxes, but 
they are substituted by the area iffi under the interpolating smooth function. In 
consequence, a curvature correction is needed since the areas added and omitted 
by this integral at the top of each box — dotted areas in the inset illustrating the 
case of one of these — are of different size. 
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Figure 1. The summation of C^ a ' accumulates the ordinate val- 
ues of a set of discrete points (bold dots). The integral under the 
smooth curve stretching from x — N + i to infinity is given by © . 



Proposition 1. 

(3) 

r(«) _ 
iV — 



dx 



=n+i/2 aTogx(logloga;) Q (a - 1) [loglog(JV + -)] ' 
Proof. This follows from the substitution 



a = 2,3, 



(4) 

via [51 2.721.2] 

(«) _ 



logx = -; x 

z 



= e 1/z ; dx = 



dz i 



fa. 
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dz 



l/log(W+l/2) 
z=0 Z(l0g^) c 



(1 — a) log" 1 z 



dz 
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1/ log(AT+l/2) 



□ 



Remark 1. The Cauchy-Maclaurin integral test plus the finiteness of this integral 
proof the convergence of 

Remark 2. f5|) estimates that an absolute accuracy o/AC' 2 ' ps 1 is reached after 
direct summation of N pa e e 1 w 15 terms, an accuracy of AC^ 2 -* ps 0.1 orafo/ a/ter 
TV fa 9 x 10 9565 terms. 
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The half-infinite interval is cut into abscissa sections of unit width centered at 
integer k, 



(G) 



(a) 
N 



E 4 Q) - 

k=N+l 



The curvature correction is estimated through a Taylor series expansion around 
the mid-point x — k in each unit interval, 



(7) A[ a) 



(8) 



dx 



k+l/2 

fc-i/2 xlogx(loglogx) a 



fe+l/2 

fc-1/2 
d 2 



f 



(x-fc) 1 



Hogfc(loglogfc) a 1! 

1 (x-fc) 2 (x-fc) 3 



dx 2 xlog x(log log x) Q 2! 



3! 



The integrals over terms with odd powers of [x — k) do not contribute due to the 
symmetry of the limits; only terms oc (x — k) 2s remain, 



(Q) A {a) = - 

V ; k fclogfc(loglogfc) 



Inserting the elementary 2.01.1] 

l.fc+1/2 



-+E 



i 



fe-i/ 



{x - kf 

2 dx 2s xlog x(log log x) a ,. (2s)! 



-dx. 



(10) 

turns this into [H Thrm. 1] 
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(x - k) 2s dx = , 



(11) A 



(a) _ 
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y7+E 



fclogfc(loglogfc) Q J 4 s (2s + 1)! dx 2s xlogx(loglogx)" 



\x—k 



The first ./V terms of ^ are summed as they stand, and the terms from TV + 1 on 
are rephrased according to (fTTj) , 
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fc=3 
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fc=JV+ 



fclogfc(loglogfc) a 



w 1 
= V - 

^fclogfcOoglogfc)- k=N+1 

to yield 

Algorithm 1. (Romberg) 
(13) 
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E K'-E 



i2» 



^ 4 s (2s + 1)! dx 2s xlogx(loglogx) Q 



logfc(log logfc)" 
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The derivatives d 2s /dx 2s share a common format, which suggests the notational 
shortcut 



4 



RICHARD J. MATHAR 



Definition 2. (Atoms of Derivatives) 
(14) 9(n,l,L) 



x n (log x) 1 (log logx) L 
The derivatives for s = 1 or s = 2 then read 
d 2 

(15) 1, = 2.9(3, 1, a) + 3. 9 (3, 2, a) + 3a 5 (3, 2, 1 + a) + 2<?(3, 3, a) 



+3ag(3, 3, 1 + a) + a(l + a)g(3, 3,1 + a), 



and 



(16) ^Lg(l, 1, a) = 24 5 (5, 1, a) + 50 5 (5, 2, a) + 70 ff (5, 3, a) + 60«?(5, 4, a) 
+24#(5, 5, a) + 50ag(5, 2, 1 + a) + 105ag(5, 3, 1 + a) + 35a(l + a)g(5, 3, 2 + a) 
+110a.g(5, 4, 1 + a) + 60a(l + a). 9 (5, 4, 2 + a) + 10a(2 + 3a + a 2 )g(5, 4, 3 + a) 
+50a 5 (5, 5, 1 + a) + 35a(l + a)g(5, 5, 2 + a) + 10a(2 + 3a + a 2 ) 5 (5, 5, 3 + a) 
+a(6 + 11a + 6a 2 + a 3 )#(5, 5,4 + a), 

for example. 

All terms are positive — so all these curvature corrections reduce the integral 
estimator jj^ in (fT3"]) . as expected from Fig. Q] for convex series of points. The 
terms of order s have a factor x 2s+1 in their denominator — they converge quicker 
than the original series. 

The properties of (|13[) for different switch-over values N are illustrated in Table 
Q] Its column is the value of (fT3"|) after replacing the infinite upper limit in 
Sfc by fc, and by including the corrections of orders s = 1, 2 and 3. The columns 

headed s = 1 or s = 2 show [4 s (2s + l)!] -1 J2k=N+i d 2s /dx 2s g(l, 1, a) to give an 
impression of the cumulative magnitude of the curvature corrections. 

Remark 3. The numbers in the column s — 2 appear to be constant as a function 
of k at the precision shown, because the contributions from the 4th derivatives are 
oc 1/fc 5 — as argued above — , so their partial sums have already converged at k = 800. 



2.2. Euler-Maclaurin. Recursive use of this technique for the sums originating 
from the curvature approximations generates an Euler-Maclaurin formula [7] 

Algorithm 2. (Centered Euler-Maclaurin) 

N oo rfts— 1 

( 17 ) cia) = E H^wg + A ^ a) + E ^w^r^. w 



Here, 



( 18 ) 2 2s -! E n V 4+2/ - I )! 

7r(s) = [l m i,2 m 2,..., s ". a ] j v J > 

is a sum over all ordered partitions (compositions) of s that accumulate the different 
paths of insertions of curvatures. We first note that the factor l/2 2s on the left 
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Table 1. Convergence of (fT31 at a = 2. 
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hand side cancels with the product FJ ■ l/(4P m i) on the right hand side. Reverse 
application of Vella's variant of Faa di Bruno's formula [HI (2)] transforms this into 



(19) (3(s) = (2^ - 1) 



(28, 



in terms of signed Bernoulli numbers B [TJ Ch. 23] with generating function [TJ 
4.5.65] 



(20) 1 - zcoshz = 2^2 (3(s) 



z 2s . 



With this approach, equations l|15 p -(|16 p become redundant and are replaced by 
derivatives of odd order in (|T7|) — to be evaluated at x = N + |. The dominant 
orders are 

(21) j-g{l, 1, a) = - 5 (2, 1, a) - g(2, 2, a) - ag(2, 2, 1 + a): 
d 3 

(22) ^3(1, 1, a) = -6<?(4, 1, a) - Ug(4, 2, a) - lla<?(4, 2, 1 + a) - 12g(4, 3, a) 
-18ag(4, 3, 1 + a) - 6a(l + a)g(4, 3, 2 + a) - 6c/(4, 4, a) - lla 5 (4, 4, 1 + a) 

-6a(l + a)ff(4, 4, 2 + a) - a (2 + 3a + a 2 )g(4, 4, 3 + a). 

Numerical examples are shown in Table [2] as a function of an upper limit s 
introduced to the s-sum in (|17[) . 



3. Summary 

The series limits of © are w 38.40676809282179, w 372.80449187938288, 
and C (4) w 3898.68733845534376. 
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Table 2. Convergence of (fTT)) . The rows s = 5 include derivatives 
up to (d 9 /dx 9 )g(l, 1, a). 
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Appendix A. The series ^2 k l/{k\og a k) 
The simpler series 

oo 1 

<23) 

can be treated by the same approach. This replaces ([3]) by 

(9 a\ j( a ) = f°° dx = I = o Q 

1 ' N L=N + i/2x(logx)" (a-l)pog(JV+l/2)]«-i' ° 

and CCH) by 

\ j V j OO 

(25) flW = £ ujZfcw + J w - E 4 «(2s + 1M ^ 



fc=2 



fe(lo g fc)« ^ z^4 s (2 s + l)! fc= ^ ' di 2 «ar(loga:) a |fc ' 
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The relevant new set of derivatives of odd order starts with 

d 1 a + log x 



dxx(\ogx) a x 2 (\ogx) 1+a ' 
d 3 1 a(2 + 3a + a 2 ) + 6alogiE + (11a + 6 log x) (log a;) 2 



(26) 
(27) 

dx 3 x(\ogx) a a; 4 (loga;) 3+Q 

The results are (rounded) w 2.10974280123689, w 2.06588653888414, 

w 2.55911974298673. The value of is known 0[1[3]. 
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